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Abstract 


A computational procedure is presented for the solution of frictional contact 
problems for aircraft tires. A Space Shuttle nose-gear tire is modeled using a two- 
dimensional laminated anisotropic shell theory which includes the effects of varia- 
tions in material and geometric parameters , transverse-shear deformation , and 
geometric nonlinearities. Contact conditions are incorporated into the formulation by 
using a perturbed Lagrangian approach with the fundamental unknowns consisting of 
the stress resultants , the generalized displacements , and the Lagrange multipliers 
associated with both contact and friction conditions . The contact- friction algorithm is 
based on a modified Coulomb friction law. A modified two-field , mixed-variational 
principle is used to obtain elemental arrays. This modification consists of augmenting 
the functional of that principle by two terms : the Lagrange multiplier vector associ- 
ated with normal and tangential node contact-load intensities and a regularization 
term that is quadratic in the Lagrange multiplier vector. These capabilities and com- 
putational features are incorporated into an in-house computer code. Experimental 
measurements were taken to define the response of the Space Shuttle nose- gear tire to 
inflation-pressure loads and to combined inflation-pressure loads and static normal 
loads against a rigid flat plate. These experimental results describe the meridional 
growth of the tire cross section caused by inflation loading , the static load-deflection 
characteristics of the tire, the geometry of the tire footprint under static loading con- 
ditions, and the normal and tangential load-intensity distributions in the tire footprint 
for the various static vertical-loading conditions. Numerical results were obtained for 
the Space Shuttle nose-gear tire subjected to inflation-pressure loads and combined 
inflation- pressure and contact loads against a rigid flat plate. The experimental mea- 
surements and the numerical results are compared. 


Introduction 

Contact-friction problems are inherently nonlinear 
and path dependent. Nonlinearity occurs partly because 
both the contact area and the contact-load intensities are 
not known beforehand and vary during the loading 
history. Path dependency is a result of the nonconserva- 
tive (irreversible dissipative) character of the frictional 
forces. 

A review of static contact problems presented in ref- 
erence 1, which includes a bibliography of approxi- 
mately 700 papers, points out that contact problems are 
important to thermomechanical stress analyses, fracture 
mechanics, mechanical problems involving elastic foun- 
dations, the mechanics of joints, geomechanics, and tires. 

Contact problems occupy a position of special 
importance in aircraft tire mechanics because the contact 
zone is where the forces are generated to support, guide, 
and maneuver the airplane. Distributions of contact loads 
and frictional forces define the moments and shears that 
are applied to the landing gear system (ref. 2). Under 
rolling conditions, the distribution of sliding velocities 
within the tire footprint combined with the frictional 
forces developed by the tire defines the rate of energy 
dissipation associated with the loading conditions and 
provides a measure of tire wear (refs. 3 and 4). In the 
case of the Space Shuttle orbiter, this wear mechanism is 


strong enough to cause tire failures during individual 
landing operations (refs. 5 and 6). Therefore, an under- 
standing of these tire friction forces and the resulting slip 
velocities is critical to the design of aircraft tires for the 
next generation of high-performance aircraft, such as the 
National Aero-Space Plane and the High-Speed Civil 
Transport. 

Modeling contact phenomena in the tire footprint is a 
formidable task partly because of difficulty of modeling 
tire response. Distribution of tractions and the footprint 
geometry are both functions of the normal, frictional, and 
inflation tire loads. Moreover, the complex mechanisms 
of dynamic friction, which allow the tire to develop the 
necessary steering and braking forces for aircraft control 
during ground operations, are not fully understood 
(ref. 7). The tire analyst thus is forced to choose among 
several friction theories. When the tire contact problem 
includes frictional effects, the solution becomes path 
dependent and a unique solution is not guaranteed. 

The aircraft tire is a composite structure of rubber 
and textile constituents that exhibit anisotropic and non- 
homogeneous material properties. Normal tire operating 
conditions create loads that can produce large deforma- 
tions. Elevated operating temperatures from the com- 
bined effects of material hysteresis and frictional heating 
can cause variations in the material characteristics of the 



tire constituents (refs. 8-10). The laminated carcass of 
the aircraft tire is thick enough to allow significant 
transverse-shear deformations. 

These facts and attendant difficulties emphasize the 
need to develop modeling strategies and analysis 
methods that include efficient, powerful and 
economical contact algorithms. Intense research has 
recently focused on nonlinear analyses of static and 
dynamic problems involving contact. Novel techniques 
that have emerged from these efforts include semi- 
analytic finite-element models for nonlinear analysis of 
shells of revolution (refs. 11 and 12), reduced methods 
(refs. 13 and 14), and operator splitting techniques 
(refs. 15-17). References 14, 17, and 18 summarize 
applications of these new tire modeling techniques. 

Objectives and Scope 

Langley tire modeling research concentrates on 
developing an accurate and efficient strategy for predict- 
ing aircraft tire responses to a variety of loading condi- 
tions. This research focuses on developing tire contact 
modeling techniques, and the specific objectives of this 
research are (1) to develop a contact algorithm with fric- 
tion effects included to predict tire response to combined 
inflation-pressure and static vertical-loading conditions, 
(2) to demonstrate the capabilities of this algorithm 
through numerical studies, and (3) to validate these 
numerical results with experimental data. Distribution of 
normal and frictional forces in the tire contact zone (or 
footprint area) is of particular interest. 

The contact algorithm is incorporated into a 
mixed-formulation, two-field, two-dimensional finite- 
element model based on the moderate-rotation 
Sanders-Budiansky shell theory, including the effects of 
transverse- shear deformations, laminated anisotropic 
material response, and nonhomogeneous shell character- 
istics (refs. 19 and 20). A perturbed Lagrangian formula- 
tion (refs. 21 and 22) is the basis for this contact 
algorithm. The Lagrangian formulation uses the pre- 
conditioned conjugate gradient (PCG) iteration proce- 
dure (refs. 23-25) to determine contact area, distribution 
of normal-force intensities, and allocation of friction- 
force intensities. A modified version of the Coulomb 
friction law is incorporated into the contact algorithm in 
which the friction coefficient at the onset of sliding is 
different from that during sliding. This algorithm also 
monitors the energy dissipated within the sliding portion 
of the contact zone. In this investigation it will be 
assumed that the tire is loaded on a surface that is much 
stiffer than the tire, thus the surface will be treated as 
rigid. Hence, the static tire contact problem will be 
treated as a unilateral contact problem. Reference 26 
summarizes the characteristics of this algorithm. 


Numerical studies presented for an inflated Space 
Shuttle nose-gear tire under static load on a flat surface 
demonstrate the capabilities of the analysis techniques. 
These analyses incorporate both friction and frictionless 
contact. Detailed studies are made of the effects of tire 
tread pattern on the contact-force intensities, the influ- 
ence of friction coefficient variations on the distribution 
of tire contact-force intensities, the convergence charac- 
teristics of the contact algorithm, and the history of 
energy dissipation in the static footprint. 

Experimental measurements were carried out on the 
Space Shuttle orbiter nose-gear tire to define its response 
to combined inflation-pressure and static vertical-loading 
conditions. This report discusses experimental proce- 
dures used to define the tire structural response to load- 
ing conditions and to measure the footprint-force 
intensities and empirical procedures used to define the 
geometry and construction details of the tire for model- 
ing purposes. Finally, the analytical results are compared 
with the experimental measurements. 

This report describes numerical studies, experimen- 
tal measurements, and comparisons between analytical 
results and experimental measurements. Reference 27 
describes development of this contact algorithm. 

Nomenclature 

Cjj y tire stiffness coefficients (/, j = 1 , 2, 6) 

d nylon cord diameter 

E |, E 2 Young’s moduli 

E , G nylon cord moduli 

E , G rubber moduli 

r r 

G 2 y G |3 , G ]2 shear moduli 

f volume fraction of nylon cord 

h total thickness of tire 

h k thickness of individual tire plies 

n normal vector to reference surface 

Po intensity of inflation pressure 

R 2 principal radii of curvature in meridional 

and circumferential directions 

r normal distance from tire axis to reference 

surface 

s , 0 shell coordinates of tire 

«, v, w displacement components of reference 

surface of tire in meridional, circumferen- 
tial, and normal directions (see fig. 1) 

w normal displacement at e = 0 = 0 (see 

fig- 7) 
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two-dimensional position vector (see 
eq. (3)) 

data vector 

Cartesian coordinate system 
extensional strains of tire reference surface 

transverse-shear strains of tire 
shear strains of tire reference surface 
Poisson’s ratio 
curvatures (see eq. (6)) 
bending strains of tire 
tension factor (see eq. (1)) 
normalized tension factor 
elongation stresses of tire 

circumferential (hoop) coordinate of tire 
(see fig. 1) 

nylon cord orientation angle (see eq. (22)) 

dimensionless coordinate along meridian 
(see fig. 4) 

rotation about normal to tire reference sur- 
face (see eq. (A9)) 

rotational components of reference surface 
of tire (see fig. 1 ) 

shear stresses of tire 


vertical-loading conditions. These sections also present 
an evaluation of the tire stiffness characteristics. 


Modeling of Tire Geometry 

A Space Shuttle nose-gear tire was cut into sections 
and used to obtain accurate measurements of the cross- 
sectional profile of an uninflated tire. A smoothed spline 
under tension was used to fit a curve through the mea- 
sured coordinates of the cross-sectional profile in a least- 
squares sense. (See refs. 28 and 29.) Because of symme- 
try, only half the cross section was modeled. A piecewise 
linear and continuous estimate of the second derivative 
c i 2 x/dz : 2 was required by the cubic spline function, and 
this was achieved by adjusting the standard deviations of 
the measured profile at the data points. The second deriv- 
ative with the effect of tension included is denoted by 


x(i') - a x(s) = 


y< : 7 

Js i+l s i 


*(s /+ i)-o y, + 1 




(i) 


for (s- <s<s i+x ) where x = [x(s), z(s)] is the position 
vector for points along the segment, y, = [jc, z]j is the 
corresponding data at point /, and a dot over a symbol 
denotes a derivative with respect to s. The chordal length 
(polygonal arc length) s i is given by the following 
equation: 


Geometric and Stiffness Characteristics of 
Space Shuttle Orbiter Nose-Gear Tire 

Tire Description 

Numerical studies conducted assess the accuracy of 
the two-dimensional shell tire model, effectiveness of the 
proposed computational procedure, and performance of 
the contact algorithm. The 32 x 8.8, type VII, 16-ply- 
rating Space Shuttle orbiter nose-gear tire was modeled 
during these studies as a two-dimensional laminated shell 
with variable thickness and variable stiffness characteris- 
tics. The outer surface of the tire was used as the refer- 
ence surface of the shell model. Figure 1 gives geometric 
characteristics of the tire. The tire carcass is constructed 
of 10 lamina of nylon and rubber with an additional rein- 
forcing ply beneath the tire tread. (See fig. 2.) The tire 
has a three-groove tread pattern. The rated load for the 
tire is 15 000 lbf at an inflation pressure of 320 psi. All 
experiments and analyses were conducted at an inflation 
pressure of 300 psi. 

The following sections describe the empirical proce- 
dures used to establish the tire geometry and to define the 
global elastic response of the tire to inflation and static 
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After solving equation (1) for x(s) and replacing x(s,) 
with y,-, the following expression is obtained: 
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Upon differentiating equation (3) and equating right- 
and left-side derivatives at s. (for i - 2, 3, I - 1), a set 
of linear algebraic equations for x(j.) is obtained. With 
the assumption of a nonperiodic spline function in which 
the slopes at and s 1 / are given, the tridiagonal differen- 
tial equation is easily solved. Once the second deriva- 
tives at point i = 1,2, / are obtained, the first and 

second derivatives at the interpolation points x(s) and 
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x(.y) are evaluated by differentiating equation (3). Refer- 
ences 29-3 1 give a detailed description of spline smooth- 
ing techniques. 

A normalized tension factor is used to eliminate a 
nonlinear behavior by setting (see ref. 32) 


a’ = 


Oisj-s^ 

l- 1 


(4) 


In practice, if this factor is less than 0.001, the resulting 
curve is approximately a cubic spline; if it is greater than 
50, the curve is nearly piecewise linear. It should be 
noted that the variable s in equations (1)~(4) is not the 
actual arc length but the chordal length approximation of 
s. Thus, the accuracy of the arc length approximation is a 
function of the number of points used in the smoothing 
procedure. 

Computed arc length coordinates x and z, and first 
and second derivatives x(^) and x(s) are used to evalu- 
ate the various geometric parameters of the tire. 

Normal vector : 


n 



Curvatures : 
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Figure 3 presents the resulting geometric characteris- 
tics of the Space Shuttle nose-gear tire. Thicknesses of 
the tire carcass at the node points of the finite-element 
model were computed along the normal vector to the tire 
reference (outer) surface by locating the points of inter- 
section of the normal vectors with the inner surface of the 
tire carcass. To facilitate these computations the tire 
inner surface was approximated by a set of third-degree 
polynomials. 


individual carcass plies was measured at the interfaces 
between the regions and these values are given in table 1 . 
A linear variation was assumed for the thickness within 
each region. Thickness of the tire tread and sidewall cov- 
ering was computed by subtracting the sum of the indi- 
vidual ply thicknesses from the total thickness of the 
carcass at each location. 

Material properties of the different plies were 
obtained with the mechanics of materials approach, 
which is widely applied to rigid composites. (See refs. 33 
and 34.) Table 2 presents elastic constants of the tire 
constituents used in this study. It was assumed that nylon 
cords of two different diameters were used in the 
construction of the tire. In region I, d = 0.022 in. for the 
bottom two plies and the tread reinforcement, and d = 
0.031 in. for all other plies. Table 3 gives cord end counts 
(epi) for individual plies at the region interfaces. A linear 
variation was assumed for epi within each region. 

Stress-strain relationships of the individual ortho- 
tropic and unidirectional layers are given by 


r 3 
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where the reduced stiffnesses c .. are given by 
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Evaluation of Stiffness Coefficients of 
Two-Dimensional Shell Model 


c 66 ~ G 12 


(14) 


The cord-rubber composite was treated as a lami- 
nated material. For the purpose of computing stiffness 
variations in the meridional direction, the tire model was 
divided into seven regions. (See fig. 2.) Thickness of the 


Elastic constants are computed by the law of mix- 
tures (see ref. 34) 

£ 1 = E c f c + E r (\-f c ) (15) 
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v 12 = U + V 1 ~fc) 

(16) 

£ r [£ c (l+2/ c ) + 2£ r (l-/ c )] 
2 ~ £ c (l-4) + 2£ r (l+0.5/ c ) 

(17) 

„ „ G r [G c + G r + (G c -G r )f c ] 

(18) 

12 " 13 “ lC c + G r -(G c -G r )f c ] 

G 23 = 0.6G 12 

(19) 

V 12 £ 2 

V,1 = 

21 £, 

(20) 


where subscripts c and r represent the quantities of the 
nylon cord and the rubber, respectively, and / is the vol- 
ume fraction of the nylon cord 

7trf 2 (epi) / 2n 

where d is cord diameter, h k is layer thickness, and epi is 
cord end count (in ends per inch). 

Stress-strain relationships of the two-dimensional 
shell were obtained by first transforming the stiffness of 
each of the individual layers to the global shell coordi- 
nates s and 0 and then integrating these coefficients 
through the thickness. Table 4 gives cord orientations in 
the individual plies of each region. The following for- 
mula was used to determine , the angle (in degrees) 
measured from the s-axis to the 0-axis, at the numerical 
quadrature points: 

0 = max[(54.382 - 3.884^ - 148.96^ 2 ) 0 , 33°] (22) 

where ^ is the dimensionless coordinate along the tire 
meridian. 

The appendix gives the resulting shell constitutive 
relations. Figure 4 shows meridional variations of the 
stiffness coefficients. 

Experimental Measurements 

A substantial experimental program has been under- 
way at Langley for several years to obtain measurements 
of the Space Shuttle nose-gear tire response to various 
loading conditions. The following sections describe these 
experimental measurements and the techniques used to 
obtain them. 


Measurements of Inflated Cross-Sectional Profile 
and Vertical Load-Deflection Response 

Close-range photogrammetry measurements . Close- 
range photogrammetry techniques were used to define 
the inflated profile of the Space Shuttle orbiter nose-gear 
tire. To facilitate these measurements, 209 circular, 
reflective targets were attached to the Space Shuttle 
nose-gear tire as shown schematically in Figure 5. Tar- 
gets were aligned along 19 meridional lines of the tire 
sidewall and a video camera was used to record the target 
positions from 10 different camera locations. A stereo- 
photography triangulation technique (refs. 35 and 36) 
was used to define the location of each target in a global 
coordinate system from these video images. The root 
mean square (rms) measurement accuracies were found 
to be 1.3 mils, 2.9 mils, and 1.5 mils in the jc-, y-, and 
^-coordinate directions, respectively. 

The same triangulation techniques were used to 
assess the ability of the loading fixture to subject the tire 
and wheel to a pure vertical load. Results of this triangu- 
lation study (illustrated in table below) indicate that there 
is a significant misalignment in the loading fixture that 
must be taken into account when analyzing the experi- 
mental results. 


Fore and aft 
displacement 

Lateral 

displacement 

Vertical 

displacement 

0.09 in. 

-0.02 in. 

1.42 in. 


Figure 6 shows unloaded and loaded profiles of a 
portion of the Space Shuttle nose-gear tire sidewall. The 
intersections of the lines in figure 6 denote the locations 
of the photogrammetry targets shown in figure 5. The 
dashed lines represent the tire sidewall geometry under a 
uniform inflation-pressure loading of 300 psi. The solid 
lines in the figure represent the tire sidewall geometry 
when the inflated tire is loaded onto a rigid, flat plate. 
For this static-loaded case, the tire cross section was 
deflected approximately 1.42 in. in the direction normal 
to the surface. Three data points are missing for the 
loaded case because the flat plate that the tire was loaded 
against interfered with the field-of-view of the camera. 
Data in figure 6 indicate that the tire sidewall in the 
vicinity of contact bulges outward radially under the 
influence of the vertical load. 

Static vertical load-deflection measurements . Static 
vertical load-deflection tests conducted on the inflated 
Space Shuttle nose-gear tire provide a global measure of 
tire elastic response. For these tests the tire was slowly 
lowered onto the flat plate until a maximum vertical load 
of approximately 30000 lbf was obtained and then 
slowly unloaded until the tire lost contact with the 
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surface. During this loading process an x-y plotter was 
used to monitor the resulting tire hysteresis loop. (See 
fig. 7.) The resulting load-deflection curve indicates that 
the nonlinear response of the tire is similar to a hardening 
spring. The hysteresis loop provides a measure of the 
damping of the tire, but for this study the damping term 
was omitted. 

Tire Footprint Geometry 

Reference 37 presents data on the footprint geometry 
for a number of tires, including the Space Shuttle 
nose-gear tire. The footprint shapes were obtained by 
applying a mixture of lightweight oil and graphite pow- 
der on the tire tread area and statically loading the tire at 
preselected vertical loads onto posterboard affixed to the 
rigid flat plate. The following paragraphs provide infor- 
mation on the measured geometric characteristics of the 
Space Shuttle nose-gear tire footprints obtained over a 
range of vertical-loading conditions at an inflation pres- 
sure of 300 psi. 

Footprint length and width measurements under 
static loading conditions . Figure 8 plots the length of the 
Space Shuttle nose-gear tire footprint as a function of 
vertical load. The symbols denote the specific incremen- 
tal load applications ranging from 2000 Ibf to 30000 lbf. 
The line through the points represents a cubic least- 
squares curve fit to the data. Figure 8 also lists the equa- 
tion representing the curve fit. Footprint length grew 
from slightly over 4 in. at a load of 2000 lbf to approxi- 
mately 14 in. at 30 000 lbf. The data indicate that the rate 
of increase in footprint length decreases as the vertical 
load increases. 

Figure 9 plots the width of the Space Shuttle 
nose-gear tire footprint as a function of vertical load. 
Various incremental loads between 2000 lbf and 
30 000 lbf are denoted by the circular symbols in 
figure 9, and the cubic least-squares curve fit to the data 
is denoted by the line through the points. The equation of 
the curve fit is included in figure 9. Footprint width 
grows from slightly more than 2 in. at the 2000-lbf nor- 
mal load case to a maximum width of about 7 in. at the 
30000-lbf loading condition. Data presented in figure 9 
indicate that the rate of growth in footprint width dimin- 
ishes as the vertical load increases. 

Figure 10 plots footprint width for the Space Shuttle 
nose-gear tire as a function of footprint length. Data pre- 
sented in this fashion indicate that there is a linear rela- 
tionship between the length and width of the Space 
Shuttle nose-gear tire over the range of vertical loads 
investigated. The least-squares equation expressing this 
relationship is shown in figure 10. 


Footprint area measurements under static loading 
conditions . Figure 1 1 plots the footprint area of the 
Space Shuttle nose-gear tire as a function of the applied 
vertical load. Numerical integration techniques were 
used to acquire the area measurements for each of the tire 
footprints at various vertical-loading conditions. For 
each footprint, two measurements were taken. The first 
measurement was obtained by computing the total area 
enclosed in the tire contact zone, including the tread 
groove area(s). These data are denoted by the circular 
symbols in figure 1 1 and are referred to as the gross foot- 
print areas. The second measurement was obtained by 
computing the contact area of the individual contact 
lobes independently. These measurements, denoted by 
the square symbols in figure 1 1, are referred to as the net 
footprint areas and exclude the tread groove area(s) for 
each footprint. Both the gross and net footprint areas 
increase in a linear fashion as the vertical load is 
increased over the range of loads investigated. The linear 
least-squares curve fit to each data set, along with its 
equation, is presented in figure 1 1 . 

Figure 12 plots the ratio of net to gross footprint 
areas for the Space Shuttle nose-gear tire as a function of 
vertical load. The data in figure 1 2 indicate that the ratio 
of net to gross footprint areas is nearly constant over the 
range of vertical-loading conditions tested. This ratio is 
approximately 0.8 for the 2000-lbf load case and about 
0.81 for the 30000-lbf load case. 

Tire Footprint Load Intensities 

Tire footprint-force transducer. Reference 37 de- 
scribes the tire footprint-force transducer used to obtain 
triaxial measurements in the contact region of the stati- 
cally loaded Space Shuttle nose-gear tire. The force 
transducer includes 10 beam-columns spaced 1.36 in. 
from center to center with a 0.25 in 2 surface area. Each 
beam-column was instrumented to measure force compo- 
nents in the normal and the two tangential directions. The 
two tangential forces represent the components of the 
friction forces in the footprint. 

The beam-columns measure total force within the 
beam-column contact surface. It should be noted that the 
instrumentation could not distinguish between distrib- 
uted loads and concentrated loads at discrete points. For 
the purpose of this investigation, the forces were 
assumed to be distributed load intensities over the net 
contact area of the individual beam-columns. 

Tire footprint areas and load intensity measure- 
ment procedures . Figure 1 3 shows schematically the tire 
footprint areas and the map of tire footprint-force trans- 
ducer locations used to obtain the data with the Space 
Shuttle nose-gear tire for three distinct vertical-loading 
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cases: 2000 lbf, 1 5 000 lbf, and 30 000 lbf. The sequence 
of testing places the tire footprint-force transducer on top 
of the rigid flat plate with the row of beam-columns 
aligned with the meridional coordinate of the tire. The 
transducer was initially positioned at approximately the 
center of the footprint. The load was then applied incre- 
mentally up to the maximum load of 30000 lbf, and data 
were recorded for each vertical-loading increment. For 
succeeding load applications, the tire footprint-force 
transducer was moved either forward or aft at half-inch 
increments in the footprint so that the entire footprint 
region was eventually covered. 

In figure 13 the map of the tire footprint-force trans- 
ducer locations is superimposed on the tire footprint. A 
trial-and-error iterative process was used to make small 
adjustments in positioning the transducer map on the tire 
footprint shapes. The final position of the transducer map 
shown in figure 13 is consistent with the pattern of non- 
zero force measurements for the various vertical-load 
cases reported in reference 37. The footprint shapes and 
the transducer outlines are drawn to scale in figure 13. It 
should be noted, however, that there is some uncertainty 
in the definition of the experimental footprint areas out- 
lines from forces such as scrubbing (or squirm) in the 
footprint areas during the application of the static vertical 
loads. 

Figure 13 also identifies the centers of contact for 
each beam-column and the percentage of actual contact 
for each beam-column that experienced less than full 
contact, i.e., loss of contact at the edge of the footprint or 
in the vicinity of tread grooves. This information was 
used to convert the transducer-force measurements into 
load intensity by dividing each normal- or tangential- 
force measurement by the appropriate beam-column con- 
tact area. The information in figure 13 also defines the 
contact boundary along the footprint periphery and along 
the edges of the tread grooves. Thus, each tire tread con- 
tact area (the ribs between the tread grooves) was treated 
as an independent zone or lobe of contact. 

Symmetries exhibited by Space Shuttle nose-gear 
tire footprint. Figure 14 shows the symmetry exhibited 
in the Space Shuttle nose-gear tire contact area. For each 
vertical-loading condition, the tire footprint is shown two 
ways. The first plot shows the mirror image of the left 
half of the tire footprint (denoted by the dashed outline) 
superimposed over the right half of the footprint. The 
second plot shows the image of the left half of the tire 
footprint (denoted by the dashed outline) rotated by 180° 
about the centroid of contact superimposed over the right 
half of the footprint. The centroids of the contact areas 
are denoted by the plus signs in figure 13. The data pre- 
sented in figure 14 indicate that the footprint of the Space 
Shuttle nose-gear tire exhibits rotational or inversion 


symmetry rather than reflective symmetry. The misalign- 
ments of the rotated dashed footprint outline relative to 
the solid footprint outline are a function of both the 
uncertainty of the footprint outline from the scrubbing 
during the loading history and imperfections in the tire 
construction. 

Contour plots. Figures 15-17 present contour plots 
of measured tire footprint load intensities in the normal 
and two tangential directions. Figure 1 3 shows the loca- 
tion of the sensors used to obtain these load-intensity 
measurements. Caution should be exercised when inter- 
preting this contour information because the contour pat- 
terns are sensitive to the sensor locations that were used 
to obtain the measurements. Note that for the 2000-lbf 
and 15 000-lbf load cases the data used to produce the 
contour plots were derived from a single row of sensors 
for each tire footprint contact lobe. For the 30000-lbf 
load case a single row of sensors was used for the two 
inner lobes and a double row of sensors was employed 
for the outer lobes. The solid lines outlining the contour 
plots represent the experimentally-determined footprint 
areas for each load case. 

Normal footprint load intensities. Figure 15 pre- 
sents the distribution of normal load intensities in the 
form of contour plots for normal loads of 2000 lbf, 
15 000 lbf, and 30000 lbf. Rated load for the Space Shut- 
tle nose-gear tire is 15 000 lbf at an inflation pressure of 
320 psi. For the 2000-lbf normal load case (fig. 15(a)) 
peak normal load intensity is between 450 and 500 psi. 
At this loading condition only the two inner tread ribs (or 
lobes) are in contact with the surface. At a normal load of 
15 000 lbf, the tire contact zone expands to include the 
two outer tread ribs. (See fig. 15(b).) Maximum normal 
load intensity for the 15 000-lbf load case is also 450 psi 
to 500 psi, but the footprint area of the tire increases sub- 
stantially over the 2000-lbf load case footprint. 

The contour plot in figure 15(c) for the 30000-lbf 
normal load condition indicates the presence of a local- 
ized maximum load intensity along the outer edge of one 
of the outside circumferential tread grooves, and is 
shown in the bottom right comer of the tire contact zone 
in figure 15(c). This peak load intensity ranges between 
650 psi and 700 psi. In the center of the outer contact 
lobes the normal load intensities range between 500 psi 
and 550 psi, and in the center of the inner lobes the nor- 
mal load intensities range between 400 psi and 450 psi 
for the 30000-lbf normal load condition. These maxi- 
mum load intensities are equal in magnitude to approxi- 
mately 150 percent to 180 percent of the inflation 
pressure. 

Tangential load intensities in lateral direction. Fig- 
ure 16 shows the lateral friction load-intensity distribu- 
tions for the Space Shuttle nose-gear tire. For the 
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2000-lbf normal load case (fig. 16 (a)), the lateral friction 
load-intensity measurements indicate negative load 
intensities of -24 psi to -30 psi on the bottom contact 
lobe and a maximum positive lateral load intensity of 
6 psi to 12 psi on the upper contact lobe. These lateral- 
traction load intensities retard the lateral growth of the 
tire footprint. Loss of skew symmetry in the lateral fric- 
tion load-intensity distribution for the 2000-lbf load case 
is attributed to the misalignment in the loading fixture 
which allowed a slight lateral displacement of the tire 
during application of the normal load. Arrows in 
figure 16(a) indicate the directions of tire deflection 
associated with this misalignment. When the normal load 
is increased to 15 000 lbf (fig. 16(b)), the magnitude of 
the peak lateral friction load intensity is also increased to 
75 psi to 90 psi and the four contact lobes exhibit a distri- 
bution of lateral load intensities which features an alter- 
nating sequence of positive and negative load-intensity 
regions. At 30000 lbf (fig. 16(c)) peak lateral friction 
load intensities are increased to 125 psi to 150 psi and the 
alternating bands of positive and negative surface trac- 
tions are especially prominent. Loss of skew symmetry 
in the lateral friction load intensities is not as obvious for 
the larger normal load cases as it is for the 2000-lbf nor- 
mal load case. Peak lateral friction load intensities occur 
near the lateral extremities of the tire footprint for the 
load conditions investigated. 

Tangential load intensities in fore and aft 
direction. Figure 1 7 shows measured fore and aft tangen- 
tial (or drag) friction load-intensity distributions for the 
Space Shuttle nose-gear tire inflated to 300 psi. For the 
2000-lbf normal load case shown in figure 17(a), drag 
friction load-intensity measurements indicate that posi- 
tive load intensities are present in both contact lobes and 
peak positive load intensities of 10 psi to 15 psi are 
located about 1 in. left of the center of contact. Fore and 
aft extremities of the tire footprint exhibit negative drag 
friction load intensities ranging up to —15 psi. Positive 
load intensities retard forward displacement of the tire 
footprint associated with misalignment of the loading 
fixture. Arrows in figure 17(a) indicate the directions of 
tire deflection associated with this misalignment. At a 
normal load of 15 000 lbf (fig. 17(b)) drag friction load 
intensities increase to a maximum positive value of 
32 psi to 40 psi in the two inner lobes about 3 in. on 
either side of the contact center. Maximum negative drag 
friction load intensities of -8 psi to —16 psi are observed 
near the right edge of contact in the two center lobes in 
figure 17(b) and there is a sharp transition between posi- 
tive and negative forces about 4 in. right of the center of 
contact. When normal force is increased to 30000 lbf 
(fig. 17(c)) the transition between positive and negative 
drag friction load intensities is very prominent. Peak drag 
friction load intensities, both positive and negative, are 


72 psi to 84 psi. About one quarter of the tire contact area 
is subjected to negative drag friction load intensities, and 
the remaining three quarters of the tire footprint are sub- 
jected to positive drag friction load intensities. Positive 
drag friction load intensities are minimum at the center 
of contact for each of the normal load cases. (See 
figs. 17(a)— (c)0 

Numerical Results and Correlation With 
Experimental Results 

Description of Finite- Element Models 

To develop the finite-element models used in the 
analysis of the Space Shuttle nose-gear tire, the cubic 
spline approximation of the outer meridional surface of 
the tire half cross section was discretized into 75 
potential node points. (See fig. 18.) From this population 
of possible nodes, a smaller number of nodes was chosen 
to approximate the tire cross section. To model the tire 
inflation response, a single strip of 30 finite elements was 
used to approximate the complete tire cross section. This 
model employed 61 nodes to characterize the tire 
meridian, and there were a total of 480 stress-resultant 
parameters and 293 nonzero generalized-displacement 
parameters to synthesize the tire inflation response. 

Finite-element models employed to analyze the con- 
tact behavior and friction characteristics of the Space 
Shuttle nose-gear tire used 41 node points in one half of a 
meridional cross section (81 nodes for the entire cross 
section) and these nodes are denoted as the circular sym- 
bols in figure 18. Nodes associated with the circumferen- 
tial tread grooves are also highlighted in figure 1 8. In the 
meridional direction, the tread area of the tire was mod- 
eled with the highest density of nodes and the sidewall 
and bead areas were modeled with progressively fewer 
nodes. This meridional node pattern was used for each of 
the two-dimensional finite-element tire models employed 
in this investigation. The circumference of the tire was 
divided into 240 possible node points and a smaller num- 
ber of nodes was chosen from that population to con- 
struct the tire finite-element models. To refine the mesh 
in specific areas such as the contact zone, a higher den- 
sity of nodes was chosen from the population in the spe- 
cific region of interest. 

Figure 19 shows a map of elements and node loca- 
tions for one of the models used to analyze the contact 
problem. Figure 19 shows an array of elements, with 40 
elements in the meridional direction and 18 elements in 
the circumferential direction. Numbers in the left and 
right margin of the figure denote the beginning and end- 
ing element numbers in specific rows. Black dots super- 
imposed over the square grid of elements denote the 
individual nodes of the finite-element model. Several 
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individual elements are shaded and shown in an 
expanded scale to illustrate the node numbering sequence 
that was used to minimize the bandwidth for the 
finite-element models. The complicating factor here is 
that the elements in the circumferential direction are 
joined along the top and bottom edge. The numbering 
scheme that is illustrated in the example shown in figure 
19 provides a minimum bandwidth for this tire model. 
For this specific example the bandwidth is 1635. The six 
rows of elements in the middle of the array, containing 
elements 1 through 240, comprise the possible contact 
region for this model. Also shown in figure 19 is the 
location of the circumferential tread grooves of the Space 
Shuttle nose-gear tire. 

Three different models were used in the analysis of 
the Space Shuttle nose-gear tire in contact with a flat 
plate. These models, denoted as model 1, model 2, and 
model 3, are depicted in figure 20. Each model employed 
480 elements in the region outside the contact zone 
(0 < -0.2 tc, 0 > 0.27t). Model 1 includes 240 elements in 
the contact region of the tire (-0.271 < 0 < 0.2 tc) for a total 
of 720 elements. (See fig. 19.) Model 1 includes 14 076 
nonzero generalized displacement parameters, 23 040 
stress resultant parameters, and 3159 contact load- 
intensity parameters. Model 2 used a refined mesh within 
the contact region, with 480 contact elements and a total 
of 960 elements. Model 2 includes 1 8 776 nonzero gener- 
alized displacement parameters, 30 720 stress resultant 
parameters, and 6075 contact load-intensity parameters. 
Model 3 employed a more refined mesh in the contact 
zone, with 960 contact elements and a total of 1440 ele- 
ments. Model 3 employed 28 152 generalized displace- 
ment parameters, 46 080 stress resultant parameters, and 

1 1 907 contact load-intensity parameters. A single itera- 
tion for model 1 required about 12 min on a Cray 2 com- 
puter, and a single iteration for model 3 required about 

1 2 min on a Cray Y-MP computer. 

Correlation Between Analytical Results and 

Experimental Measurements 

Tire response to inflation-pressure loading. To 
assess the accuracy of the high fidelity carcass model of 
the Space Shuttle nose-gear tire, profile deformations 
produced by a uniform inflation pressure of 300 psi and 
acting normal to the interior surface were calculated 
using a strip of geometrically nonlinear shell finite 
elements described above. Figure 21 presents inflated 
and uninflated cross-sectional profiles for the Space 
Shuttle nose-gear tire. The primary effect of inflation 
pressure is to expand the tire profile in the cross-sectional 
regions I-V. (See fig. 2.) The predicted inflated profile is 
in excellent agreement with the measured profiles from 
the photogrammetry studies. Photogrammetry data from 


a number of tire cross sections (denoted by the different 
symbols) are compared with the model predictions in fig- 
ure 21. These results clearly illustrate the axisymmetric 
characteristics of the tire response to inflation loads. Ref- 
erence 38 presents additional information on inflation 
pressure results. 

Tire load-deflection and load-contact-area charac- 
teristics, One global measure of tire response to static 
contact loading is the normal load-deflection curve. Fig- 
ure 22 presents measured and calculated normal load- 
deflection curves for the Space Shuttle nose-gear tire 
inflated to 300 psi and loaded against a flat plate. The 
measured load-deflection curve is denoted by square 
symbols in figure 22 and a quadratic curve is used to fair 
the measurements in a least-squares fashion. The calcu- 
lated load-deflection response of the tire was obtained by 
taking the product of the predicted normal load intensi- 
ties and the calculated contact area at each contact node 
and summing them over all contact nodes. Circular sym- 
bols denote the predicted load-deflection response based 
on the Newton -Cotes quadrature calculation of the 
contact areas associated with each contact node. Triangu- 
lar symbols denote the predicted load-deflection 
response based on the contact area algorithm described in 
reference 27. Both calculation procedures perform well 
in matching the measured normal tire load-deflection 
response, but the calculation procedure based on the con- 
tact area algorithm appears to be more accurate over the 
load range analyzed. 

A second global measure of tire response to static 
contact loading is the normal-load-contact-area behavior. 
Figures 23 and 24 present measured and calculated 
normal-load-contact-area characteristics of the Space 
Shuttle nose-gear tire subjected to static contact loading 
against a flat plate. In figure 23 measured and predicted 
net footprint areas are plotted as a function of the applied 
normal load. Net footprint area is defined as the sum of 
the contact areas of the individual tire tread ribs or lobes 
(see fig. 15, for example) and excludes the areas associ- 
ated with the circumferential grooves. In figure 24 the 
gross footprint areas are plotted as a function of the 
applied normal load. Gross footprint area is defined as 
the total contact area of the tire, and the groove areas are 
included in this parameter. For a smooth tread tire, the 
net and gross footprint areas are identical. 

Figures 23 and 24 denote experimental data by 
square symbols. Predicted contact areas based upon the 
Newton-Cotes quadrature calculation are denoted by cir- 
cular symbols and the dashed linear curve representing a 
least-squares fairing of the analytical results. Predicted 
contact areas based upon the contact- area algorithm cal- 
culations are denoted by triangular symbols and the solid 
linear curve least-squares fairing of the analytical results. 
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Results presented in figures 23 and 24 indicate that 
contact-area predictions based upon the contact-area 
algorithm are much more accurate than those based upon 
the Newton-Cotes quadrature calculation. Predicted net 
footprint areas are more accurate than those of the gross 
footprint areas regardless of the calculation method used. 
Results presented in figures 23 and 24 also indicate that 
additional fine-tuning of the contact-area algorithm is 
needed for improved contact-area predictions, especially 
at the higher normal load conditions. 

Figure 25 presents the method of obtaining a predic- 
tion of the tire footprint shape from the analytical results. 
Figure 25 shows a map of contact nodes, denoted by 
square symbols, in and around the footprint area pro- 
jected onto the contact surface. Dark symbols represent 
nodes that are in contact with the flat surface and light 
symbols represent nodes that are not in contact. A solid 
line denotes the measured footprint outline obtained from 
the lightweight oil and graphite powder footprint stamps 
described previously in reference 37. A dashed line rep- 
resents an elliptical curve approximation of the footprint 
shape based on the distribution of contact nodes. For the 
2000-lbf normal load case shown in figure 25(a), mea- 
sured and calculated footprint widths are in close agree- 
ment, but predicted footprint length is less than the 
measured length. These measured and/or calculated foot- 
print shapes are included in figures 26-31 as an aid in 
interpreting the footprint load-intensity data. For the 
15 000-lbf normal load case shown in figure 25(b), mea- 
sured and calculated footprint shapes are in close agree- 
ment. It should be noted that the footprint outline looks 
distorted in figure 25(b) due to a difference in scaling 
factors between the two coordinates. 

Tire contact load-intensity distributions. Fig- 
ures 26-28 present comparisons of measured and pre- 
dicted tire footprint load-intensity distributions. Data 
presented in the form of contour plots and measured 
load-intensity distributions shown in figures 15-17 are 
reproduced here to aid in the comparison. Boundary con- 
ditions used to produce the analytical results were chosen 
to simulate the misalignment in the loading fixture 
described earlier and arrows in figures 26-28 denote the 
direction, but not the magnitude, of tangential displace- 
ments that occurred during the normal loading sequence 
as a result of the misalignment. Part (a) of each figure 
presents results from the 2000-lbf normal load case, 
and part (b) of each figure presents results from the 
1 5 000-lbf load case. All analytical results shown in fig- 
ures 26-28 were obtained using model 1 . 

For the 2000-lbf normal load case the analysis 
predicts a ridge of high normal load intensities along 
both edges of the center circumferential tread groove. 
(See fig. 26(a).) The magnitude of normal load intensities 


along these ridges is between 900 psi and 1000 psi. 
Experimental measurements of the normal load-intensity 
distribution cannot confirm the existence of these high 
load-intensity ridges since no load transducers were 
located in the critical region along the edge of the center 
groove for this normal load condition. For the remainder 
of the footprint area the predicted magnitudes of normal 
load intensities are slightly less than the measured inten- 
sities, and the predicted distribution of normal load inten- 
sities is generally in fair agreement with experimental 
measurements. These results are consistent with analyti- 
cal results for a tire tread presented in reference 39. 

For the 15 000-lbf normal load case the analysis pre- 
dicts ridges of high load intensities along each edge of 
the three circumferential tread grooves. Magnitudes 
associated with these ridges range between 500 psi and 
800 psi, and are generally lower than the peak load inten- 
sities predicted for the 2000-lbf normal load case. Pre- 
dicted and measured normal load-intensity distributions 
are in good agreement. 

Figure 27 shows comparison of measured and pre- 
dicted lateral friction load-intensity distributions. For the 
2000-lbf normal load case (fig. 27(a)), both measured 
and predicted lateral friction load-intensity distributions 
exhibit loss in skew symmetry across the width of the 
footprint due to misalignment in the loading fixture. 
Measured and calculated lateral friction load-intensity 
distributions are in good agreement, and the magnitudes 
of the measured and predicted lateral friction load inten- 
sities are in fair agreement. 

For the 15000-lbf normal load case (fig. 27(b)), both 
predicted and measured lateral friction load intensities 
are shown to reach their maximum magnitudes in the lat- 
eral extremities of the tire footprint. This trend is more 
apparent for predicted results than for measured values. 
Both predicted lateral friction load intensities and 
measured results exhibit bands of alternating positive and 
negative friction values across the width of the tire 
footprint. 

Figure 28 shows comparison of measured and pre- 
dicted drag friction load-intensity distributions. For the 
2000-lbf normal load case (fig. 28(a)) predicted drag 
friction load intensities are higher than the measured val- 
ues. Experimental measurements and analytical predic- 
tions both indicate local maximums in the drag friction 
load intensities in the fore and aft sections of the tire 
footprint. However, the small negative drag friction load 
intensities seen in the experimental results are not present 
in the analytical predictions. Overall, the comparison 
between measured and predicted drag friction load- 
intensity distributions indicate that analytical predictions 
give a reasonable approximation of the drag friction dis- 
tribution of the tire for static loading conditions, but 
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additional model fine-tuning efforts are needed to bring 
analytical and experimental drag friction load-intensity 
results into closer agreement. 

For the 15 000-lbf normal load case (fig. 28(b)), both 
predicted and measured drag friction load-intensity dis- 
tributions are shown to exhibit local positive and nega- 
tive maximums in the fore and aft sections of the tire 
footprint and to exhibit minimum magnitudes near the 
center of the footprint. Predicted drag friction load- 
intensity magnitudes are generally higher than measured 
values. The analysis predicts the formation of larger 
pockets of negative drag friction load intensities in the 
leading edge of the footprint (the right edge of the foot- 
print in fig. 28(b)) than the measurements indicate, but 
this trend is consistent with the experimental measure- 
ments obtained for the 30 000-lbf normal load case 
shown in figure 17(c). 

Correlation between the analytical predictions and 
experimental measurements of tire footprint area, normal 
load-intensity distributions, and tangential load-intensity 
distributions appears to be stronger for the 15 000-lbf 
normal load case than for the 2000-lbf load case. 

Influence of Model Refinement on Predicted 

Contact Load-Intensity Distributions 

The contour plots in figure 29 show predicted con- 
tact load-intensity distributions for the Space Shuttle 
nose-gear tire subjected to an inflation pressure of 
300 psi and unsymmetric static loading conditions asso- 
ciated with the experimental loading fixture misalign- 
ment. Load-intensity distribution predictions were made 
with model 1 and model 3 (figs. 19 and 20) for the 
2000 lbf normal load case. Tangential displacement 
directions associated with the loading fixture misalign- 
ment are denoted by arrows in figure 29. Calculated con- 
tact area for each model is denoted by a dashed outline in 
each contour plot and there is no significant difference 
between the two calculated footprint shapes. Figure 29(a) 
shows calculated normal load-intensity distributions. The 
more refined model 3 and model 1 each predicted the 
existence of a ridge of high load intensities along both 
edges of the center tread groove for this loading case, but 
the predicted magnitudes of these load intensities from 
model 3 were about 200 psi lower than those predicted 
by model 1. This result is consistent with results pre- 
sented in reference 38 for frictionless contact. Overall, 
the predicted normal load-intensity distribution from 
model 3 is smoother than the predicted distribution from 
model 1. 

Figure 29(b) presents calculated lateral friction load- 
intensity distributions for model 1 and model 3. Gener- 
ally, the predicted lateral friction load-intensity magni- 
tudes from model 3 are larger than those from model 1 


for the loading case considered. Both models predict a 
loss in lateral friction load-intensity skew symmetry due 
to the unsymmetric loading condition and this lack 
of symmetry is more pronounced for the model 1 
predictions. 

Figure 29(c) presents calculated drag friction load- 
intensity distributions for the models. Predicted drag fric- 
tion load intensities from model 3 are consistently higher 
than those from model 1, and the model 3 drag friction 
load-intensity distribution is more uniform over the tire 
contact region than the predicted drag friction load- 
intensity distribution from model 1 . 

In summary, predicted contact load-intensity distri- 
butions for model 1 and model 3 are generally in close 
agreement. Model 3 softens the predicted peak normal 
load intensities along the edges of the center tread 
groove, but predicts larger friction load-intensity magni- 
tudes than model 1. Predicted normal load-intensity 
distribution from model 3 is in closer agreement with the 
experimental measurements (see figs. 15-17 and 
figs. 26-28 for examples) than that from model 1 , but 
predicted lateral friction and drag friction load-intensity 
distributions from model 1 are closer to the experimental 
measurements than those from model 3. Both models 
predict a shorter tire footprint length than the experimen- 
tal measurement for the 2000-lbf normal load case. 

Effect of Tread Grooves and Friction on 

Predicted Contact Load-Intensity Distributions 

Figure 30 presents a comparison of predicted contact 
load-intensity distributions for a smooth-tread Space 
Shuttle nose-gear tire and a circumferentially-grooved 
tread Space Shuttle nose-gear tire subjected to a normal 
load of 2000 lbf. Analysis of both tread configurations 
was conducted with model 1, and for these results, tire 
normal loading was assumed to be symmetric with no 
tangential displacements associated with load fixture 
misalignments. Figure 30(a) shows normal load-intensity 
distributions. Smooth-tread tire results indicate peak nor- 
mal load intensities between 300 psi and 400 psi near the 
center of contact. When the smooth tread is replaced with 
a circumferentially-grooved tread, two ridges of high 
normal load intensities are observed along the edges of 
the center groove. Magnitude of normal load intensities 
along these ridges is between 900 psi and 1000 psi. It 
should be noted that normal load-intensity distribution 
for the symmetrically- loaded grooved-tread tire is nearly 
identical to the predicted normal load-intensity distribu- 
tion shown in figure 26(a) for the unsymmetrically 
loaded tire. The presence of the circumferential tread 
grooves also affects the predicted footprint shape. The 
grooved-tread footprint is slightly longer and thinner 
than the smooth-tread footprint. 


11 



Figure 30(b) shows comparison of the lateral friction 
load-intensity distributions for the grooved- and smooth- 
tread tire and figure 30(c) shows the drag friction load- 
intensity distribution comparison. Lateral friction load 
intensities for both tread configurations exhibit skew 
symmetric distributions across the width of the tire foot- 
print and the presence of the center circumferential 
groove does not have a significant effect on the lateral 
friction load-intensity distribution. Insensitivity of the 
lateral friction load-intensity distribution to the presence 
of a center groove is attributed to the natural distribution 
of lateral friction load intensities which concentrate these 
friction forces in the lateral extremities of the footprint. 
Predicted drag friction load intensities for this symmetric 
loading condition are uniformly small over the entire 
footprint and should be compared with the predicted drag 
friction load intensities for the unsymmetric loading 
case. (See fig. 28(a).) 

In summary, a comparison of the predicted load- 
intensity distributions presented in figure 30 with those 
presented in figures 26-28 indicates that predicted nor- 
mal load-intensity distributions are strongly influenced 
by the presence of tread grooves and generally un- 
affected by unsymmetric loading conditions associated 
with the loading fixture misalignments. Conversely, lat- 
eral friction and drag friction load-intensity distributions 
were unaffected by the presence of tread grooves for the 
loading conditions examined, but were very sensitive to 
tangential displacements associated with the misalign- 
ment in the loading fixture. 

Figure 31 presents a comparison of predicted normal 
load intensities for frictional and frictionless contact of a 
grooved-tread Space Shuttle nose-gear tire. Both analy- 
ses were conducted with model 1 . The comparison pre- 
sented in figure 31 indicates that the normal load- 
intensity distribution for symmetric static loading cases 
is unaffected by the presence of friction forces. This 
result implies that a reasonable modeling strategy for 
tires subjected to static loading contact conditions might 
involve operator splitting techniques to first apply the 
normal load with a frictionless contact assumption and 
then apply subsequent tangential loads with a contact- 
friction algorithm. 

Strain Energy Density Distributions 

One method of showing the regions of high strain 
due to the imposed loads on a complex structure such as 
a tire is strain energy density. Calculated variations in the 
strain energy density for the Space Shuttle nose-gear tire 
are presented in figure 32 in the form of contour plots. 
Results are shown for normal loading cases of 2000 lbf 
and 15000 lbf. Figure 32(a) presents total strain energy 
density and figure 32(b) presents transverse- shear strain 


energy density. For the load cases shown, the total strain 
energy density is primarily influenced by the inflation 
pressure load. Since the normal loads associated with 
contact are compressive, the total strain energy density is 
reduced in the region of contact. Total strain energy den- 
sity is also reduced in the lower tire sidewall near the 
bead. Transverse-shear strain energy density is maxi- 
mized in the tire sidewall near the contact zone and along 
the lower sidewall near the bead. In the upper sidewall 
and shoulder of the tire near the contact region (see 
fig. 32(b) insert) the transverse-shear strain energy repre- 
sents about 25 percent to 50 percent of the total strain 
energy associated with static contact. 

Conclusions 

A computational procedure is presented for the solu- 
tion of frictional contact problems for aircraft tires. The 
Space Shuttle nose-gear tire is modeled using a two- 
dimensional laminated anisotropic shell theory with the 
effects of variations in material and geometric parame- 
ters, transverse-shear deformation, and geometric non- 
linearities included. 

Experimental results are presented that describe the 
response of the Space Shuttle nose-gear tire to inflation- 
pressure loads and to combined inflation-pressure loads 
and normal static loads against a rigid flat plate. These 
experimental results describe the meridional growth of 
the tire cross section due to inflation loading, the static 
load-deflection characteristics of the tire, the geometry of 
the tire footprint under static loading conditions, and the 
normal and tangential load-intensity distributions in the 
tire footprint for the various static vertical -loading 
conditions. 

Numerical results are presented for the Space Shuttle 
nose-gear tire subjected to inflation-pressure loads and 
combined inflation-pressure and contact loads against a 
rigid flat plate. 

Results from this investigation lead to the following 
observations and conclusions: 

1 . The predicted inflated profile is in excellent agree- 
ment with measured profiles from the photogram- 
metry studies. Predicted normal and tangential 
contact load-intensity distributions are in fair to 
good agreement with experimental measurements, 
though additional experimental and analytical stud- 
ies are needed. Correlation between analytical pre- 
dictions and experimental measurements of tire 
footprint area, normal load-intensity distributions, 
and tangential load-intensity distributions appears 
to be stronger for the 15 000-lbf normal load case 
than for the 2000-lbf load case. 
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2. Predicted normal load-intensity distribution from 
the model with a refined grid in the tire contact area 
is in closer agreement with the experimental mea- 
surements than the distribution from the least 
refined model used in this investigation, but pre- 
dicted lateral friction and drag friction load- 
intensity distributions from that least refined model 
are closer to the experimental measurements than 
those from the refined model. 

3. Predicted normal load-intensity distributions are 
strongly influenced by the presence of tread 
grooves and generally unaffected by unsymmetric 
loading conditions associated with the loading fix- 
ture misalignments. Lateral friction and drag fric- 
tion load-intensity distributions are unaffected by 
the presence of tread grooves for the loading condi- 
tions examined, but are very sensitive to tangential 
displacements associated with the misalignment in 
the loading fixture. 


4. Normal load-intensity distributions for symmetric 
static loading cases are unaffected by the presence 
of friction forces. This result implies that a reason- 
able modeling strategy for tires subjected to static- 
loading contact conditions might involve operator 
splitting techniques to first apply the normal 
load with a frictionless contact assumption and 
then apply subsequent tangential loads with a 
contact-friction algorithm. 

5. In the upper sidewall and shoulder of the tire near 
the contact region the transverse-shear strain 
energy represents about 25 percent to 50 percent of 
the total strain energy associated with static 
contact. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
February 6, 1996 
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Appendix 


Fundamental Equations of Shell Theory Used in Present Study 

This appendix summarizes the fundamental equations of the Sanders-Budiansky type shell of revolution used in this 
study. Effects of laminated, anisotropic material response and transverse-shear deformation are included in these 
relationships. 


Strain-Displacement Relationships 
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where e s and e 0 are extensional strains in the meridional and circumferential directions, 2e i0 is the in-plane shear 
strain, K^and K 0 are bending strains in the meridional and circumferential directions, 2 k 50 is the twisting strain, 2e v3 
and 2e 03 are transverse-shear strains, d s = g- , d 0 = , and <j> is the rotation around the normal to the shell, which is 

given by 4 0 


♦ = 


— d a u + f 3 



} 


Nonlinear terms that account for moderate rotations are underlined with dashes in equations (Al) -(A3). 


(A9) 


Constitutive Relations 

The shell is assumed to be made of a laminated, anisotropic, linearly elastic material. Every point of the shell is 
assumed to possess a single plane of elastic symmetry parallel to the middle surface. The relationships between the 
stress resultants and the strain measures of the shell are given by 
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(A 10) 


where c (y , /-, and d tJ ( i,j = 1, 2, 6) are shell stiffness coefficients. Nonorthotropic (anisotropic) terms are circled and dots 
indicate zero terms. 
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Table 1 . Variation of Ply Thickness hiJh 0 
[Iiq = 0.7513 in.] 


Ply number 
(top to bottom) 

V*0 f° r re gi° n — 

I 

II 

III 

IV 

V 

VI 

VII 

i 

(tread and sidewall) 

16 

h \ = h ~Y* h k 

k= 2 

a 2 

0.0865- 

.0865 

0.0865- 

.0658 

0.0658- 

.0692 

0.0692- 

.0813 

0.0801- 

.0937 

0.0681- 

.1238 

0.0918- 

.1240 

3 

0.0865- 

.0865 

0.0865- 

.0658 

0.0658- 

.0692 

0.0692- 

.0813 

0.0801- 

.0937 

0.0681- 

.1238 

0.0918- 

.1240 

4 

0.0865- 

.0865 

0.0666- 

.0506 

0.0506- 

.0532 

0.0692- 

.0813 

0.0488- 

.0571 

0.0523- 

.0950 

0.0652- 

.0880 

5 

0.0666- 

.0666 

0.0666- 

.0506 

0.0506- 

.0532 

0.0426- 

.0500 

0.0488- 

.0571 

0.0523- 

.0950 

0.0652- 

.0880 

6 

0.0666- 

.0666 

0.0666- 

.0506 

0.0506- 

.0532 

0.0426- 

.0500 

0.0488- 

.0571 

0.0523- 

.0950 

b 0.2662- 

.3594 

7 

0.0666- 

.0666 

0.0666- 

.0506 

0.0506- 

.0532 

0.0426- 

.0500 

0.0488- 

.0571 

0.0523- 

.0950 

0.0652- 

.0880 

8 

0.0666- 

.0666 

0.0666- 

.0506 

0.0506- 

.0532 

0.0426- 

.0500 

0.0488- 

.0571 

0.0523- 

.0950 

0.0652- 

.0880 

9 

0.0666- 

.0666 

0.0666- 

.0506 

0.0506- 

.0532 

0.0426- 

.0500 

0.0488- 

.0571 

0.0523- 

.0950 

0.0652- 

.0880 

10 

0.0666- 

.0666 

0.0466- 

.0354 

0.0354- 

.0373 

0.0426- 

.0500 

0.0488- 

.0571 

0.0523- 

.0950 

0.0652- 

.0880 

11 

0.0466- 

.0466 

0.0466- 

.0354 

0.0354- 

.0373 

0.0346- 

.0407 

0.0488- 

.0571 

0.0523- 

.0950 

b 0.2662- 

.3594 

12 

0.0466- 

.0466 

0.0798- 

.0798 

0.0798- 

.0798 

0.0346- 

.0407 

0.0375- 

.0439 

0.0523- 

.0950 

0.0652- 

.0880 

13 

0.0798- 

.0798 

0 

0 

0.0798- 

.0798 

0.0375- 

.0439 

0.0523- 

.0950 

0.0652- 

.0880 

14 

0 

0 

0 

0 

0.0798- 

.0798 

0.0366- 

.0666 

0.0466- 

.0629 

15 

0 

0 

0 

0 

0 

0.0366- 

.0666 

0.0466- 

.0629 

16 

0 

0 

0 

0 

0 

0.0798- 

.1464 

0.1464- 

.1597 


a Second layer of region I represents layer which has reinforcement (see fig. 2). 
Represents thickness of bead wires. 
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Table 2. Values of Elastic Constants of Tire Constituents Used in Present Study 


Tire constituent 

Young’s modulus, psi 

Shear modulus, psi 

Poisson’s ratio 

Rubber 

4.5 x 10 2 

1.51 x 10 2 

0.49 

Nylon cord 

3.5 x 10 5 

7.00 x 10 2 

.66 

Bead a 

2.9 x 10 7 

1 . 1 0 x 1 0 7 

.30 


a Since deformations are small in bead area, it is reasonable to assume that bead wires are isotropic. 


Table 3. Variation of Nylon Cord End Counts in Different Plies Along Meridian 



Cord end count, ends per inch, for region — 

Ply number 
(top to bottom) 

i 

II 

III 


V 

VI 

VII 

1 

(tread and sidewall) 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

2 

16-16 

18-14 

14-14 

14-14 

14-14 

14-14 

14-14 

3 

18-18 

18-14 

14-14 

14-14 

14-14 

14-14 

14-14 

4 

18-18 

21-20 

20-18 

18-16 

16-16 

16-14 

14-14 

5 

23-21 

21-20 

20-18 

18-16 

16-16 

16-14 

14-14 

6 

23-21 

21-20 

20-18 

18-16 

16-16 

16-14 

Bead 

7-9 

23-21 

21-20 

20-18 

18-16 

16-16 

16-14 

14-14 

10 

23-21 

29-26 

26-25 

18-16 

16-16 

16-14 

14-14 

11 

30-29 

29-26 

26-25 

25-24 

16-16 

16-14 

Bead 

12 

30-29 

Rubber 

Rubber 

25-24 

24-22 

16-14 

14-14 

13 

Rubber 



Rubber 

24-22 

16-14 

14-14 

14 





Rubber 

22-22 

22-22 

15 






22-22 

22-22 

16 






Rubber 

Rubber 
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Table 4. Variation of Cord Orientation of Individual Plies, 0* , Along Meridian 





0* , deg* for region- 




Ply number 
(top to bottom) 

I 

II 

III 

IV 

V 

VI 

VII 

1 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 





_ 




2 

-0-6 

0 

0 

0 

0 

0 

0 



_ 



_ 



3 

0 

-0 

-0 

-0 

-0 

-0 

-0 









4 

-0 

0 

0 

-0 

-0 

-0 

0 









5 

0 

-0 

-0 

0 

0 

0 

-0 

6 

-0 

0 

0 

-0 

0 

-0 

Bead 




__ 

_ 




7 

0 

-0 

-0 

0 

-0 

0 

-0 

8 



_ 





-0 

G 

0 

-0 

0 

-0 

0 









9 

0 

-0 

-0 

0 

-0 

0 

-0 

10 








-0 

0 + 6 

0 + 6 

-0 

0 

-0 

0 

11 

0 + 6 

-0-6 

-0-6 

0 + 6 

-0 

0 

Bead 

12 

13 

14 

15 

16 

-0-6 

Rubber 

Rubber 

Rubber 

-0-6 

Rubber 

0 + 6 

-0-6 

Rubber 

0 

-0 
0 + 6 
-0-6 
Rubber 

0 

-0 
0 + 6 
-0-6 
Rubber 


20 













Boundary conditions: 

At £ = 0.5 At £ = 0 

u = v = w = 0 Symmetry condition: 

<t> s = <l> 6 = 0 u = v = <t> s = <j) e = 0 


Figure 1 . Geometric characteristics of Space Shuttle nose-gear tire. 
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Normalized 

parameter 




Figure 3. Meridional variation of geometric parameters of two-dimensional shell model of Space Shuttle nose-gear tire. 
Reference surface chosen to be outer surface. r 0 = 15.1737 in.; k 2? o = 0.0659 in -1 ; Kj 0 = 0.1091 in" 1 . 
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Vertical 

load, 

kips 



Deflection, w, in. 


Figure 7. Static vertical load-deflection curve for Space Shuttle nose-gear tire. 



Figure 8. Footprint length as a function of vertical load for Space Shuttle nose-gear tire. Inflation pressure = 300 psi 





1 h 


I 1 1 I I 

0 10 20 30 40 x 10 3 

Vertical load, Ibf 


Figure 9. Footprint width as a function of vertical load for Space Shuttle nose-gear tire. Inflation pressure = 300 psi. 


Footprint 

width, 

in. 



Figure 10. Footprint width as a function of footprint length for Space Shuttle nose-gear tire. Inflation pressure = 300 psi. 


80 


60 


Footprint 

area, 

in^ 


40 


20 
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Left half of footprint reflected 
onto right half of footprint 


Left half of footprint rotated 
onto right half of footprint 



(a) Vertical load = 2000 lbf. 

Figure 14. Symmetries exhibited by footprint of Space Shuttle nose-gear tire. 
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Tangential displacement directions 


1 



Footprint length, in. 


(a) Normal load, 2000 lbf. 



Footprint length, in. 

(b) Normal load, 15000 lbf. 


Normal load intensity, psi 


800 

700 

600 

500 

400 

300 

200 

100 

0 



Footprint length, in. 

(c) Normal load, 30000 lbf. 

Figure 15. Color contour plot of measured footprint normal load-intensity distribution for Space Shuttle orbiter nose- 
gear tire subjected to an inflation pressure of 300 psi. 
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Tangential displacement directions 


I 



Footprint length, in. 


Lateral friction load intensity, psi 


42 

30 

18 

6 


48 

36 

24 

12 


-18 

-30 

-42 


-12 

-24 

-36 

-48 


(a) Lateral friction load intensity for a normal load of 2000 lbf. 



Footprint length, in. 


Lateral friction load intensity, psi 


105 

75 

45 

15 


-15 


-45 


-75 

-105 


120 

90 

60 

30 

0 


-30 


-60 


-90 

-120 


(b) Lateral friction load intensity for a normal load of 15000 lbf. 


Lateral friction load intensity, psi 



Footprint length, in. 


175 

125 

75 

25 


-25 

-75 

-125 

-175 


200 

150 

100 

50 

0 

-50 

-100 


-150 


-200 


(c) Lateral friction load intensity for a normal load of 30000 lbf. 


Figure 16. Color contour plot of measured footprint lateral 
gear tire subjected to an inflation pressure of 300 psi. 


load-intensity distribution for Space Shuttle orbiter nose- 
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Tangential displacement directions 


I 



Footprint length, in. 


Drag friction load intensity, psi 


35 

25 

15 

5 

-5 

-15 

-25 

-35 


40 

30 

20 

10 

0 

-10 

-20 

-30 

-40 


(a) Drag friction load intensity for a normal load of 2000 lbf. 



Footprint length, in. 


Drag friction load intensity, psi 


84 

60 

36 

12 

-12 

-36 

-60 

-84 


96 

72 
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0 

-24 
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(b) Drag friction load intensity for a normal load of 15000 lbf. 



Footprint length, in. 


Drag friction load intensity, psi 


84 

60 


36 H 


12 


-12 


-36 


-60 

-84 


96 

72 

48 

24 

0 


-24 

-48 


-72 
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(c) Drag friction load intensity for a normal load of 30000 lbf. 


Figure 17. Color contour plot of measured footprint drag load-intensity distribution for Space Shuttle orbiter nose-gear 
tire subjected to an inflation pressure of 300 psi. 
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Figure 1 8. Meridional profile of Space Shuttle nose-gear tire models denoting node point locations. 





Figure 19. Typical array of finite elements and nodes used to model Space Shuttle nose-gear tire. 
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Model 1 


Model 2 


Model 3 

i 



M9-4 finite element 


Sector 

Number of elements (M9-4) for — j 

Model 1 

Model 2 

Model 3 

-0.27: < 0 < 0.2ti 

240 (40 x 6) 

480 (40 x 12) 

960 (40 x 24) 

0 < -0.2ti, 0 > 0.2n 

480 (40 x 12) 

480 (40 x 12) 

480 (40 x 12) 

Total 

720 (40 x 18) 

960 (30 x 24) 

1440 (40 x 36) 


Figure 20. Finite-element models of Space Shuttle nose-gear tire used in present study. 


Internal inflation loading, 
300 psi 




Figure 2 1 . Inflated and uninflated profiles of Space Shuttle nose-gear tire. 
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Tire normal 
load, lbf 


— 

□ 

Experimental measurement 


A 

Contact area algorithm calculation 


O 

Newton-Cotes quadrature calculation 



0 .5 1.0 1.5 2*0 2.5 3.0 

Tire normal deflection, in. 


Figure 22. Comparison of calculated and measured normal load-deflection curve for Space Shuttle nose-gear tire. Infla- 
tion pressure = 300 psi. 
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0 10 000 20 000 30 000 40 000 

Tire normal load, lbf 


Figure 23. Comparison of calculated and measured net tire footprint areas for Space Shuttle nose-gear tire. Inflation 
pressure = 300 psi. 
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□ 

Experimental measurement 

— 

A 

Contact area algorithm calculation 

— 

o 

Newton-Cotes quadrature calculation 


Gross footprint 
area, in 2 



Figure 24. Comparison of calculated and measured gross tire footprint areas for Space Shuttle nose-gear tire. Inflation 
pressure = 300 psi. 
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Footprint 
width, in. 



Footprint length, in. 

(a) Normal load = 2000 lbf. 


Footprint 
width, in. 



-8 -6 -4 -2 0 2 4 6 8 


Footprint length, in. 


(b) Normal load = 15000 lbf. 


Figure 25. Measured and calculated footprints for Space Shuttle nose-gear tire. Inflation pressure = 300 psi; analytical 
results from model 1 . 
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Tangential displacement directions 


Experimental measurement 



Footprint length, in. 


Analytical prediction 



Footprint length, in. 


Normal load 
intensity, psi 


1000 

900 

800 

700 

600 

500 

400 

300 

200 

100 

0 


(a) Normal load = 2000 lbf. 




Footprint length, in. Footprint length, in. 

(b) Normal load = 15 000 lbf. 


Figure 26. Comparison of measured and predicted footprint normal load-intensity distribution for Space Shuttle 
nose-gear tire subjected to an inflation pressure of 300 psi. For model 1: p static = 0.6; ^dynamic = 0.51; E n = 1.0E+12; 
£, = 2.0E+03; Creiax = 0.5. 
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Tangential displacement directions 




Experimental measurement 



Footprint length, in. 


Analytical prediction 


Lateral friction 
load intensity, psi 
48 
42 

36 




Footprint length, in. 


(a) Lateral friction load intensity for a normal load of 2000 lbf. 



Footprint length, in. 


Lateral friction 
load intensity, psi 



120 

90 

60 

30 

0 

-30 

-60 

-90 

120 


Footprint length, in. 


(b) Lateral friction load intensity for a normal load of 15000 lbf. 


Figure 27. Comparison of measured and predicted footprint lateral friction load-intensity distribution for Space Shuttle 
nose-gear tire subjected to an inflation pressure of 300 psi. For model 1: p static = 0.6; Adynamic = 0.51; e„ = l.OE+12; 
= 2.0E+03; £re] a x = 
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Tangential displacement directions 


I 


Experimental measurement 



Footprint length, in. 


Drag friction 
load intensity, psi 
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(b) Drag friction load intensity for a normal load of 15000 lbf. 


Figure 28. Comparison of measured and predicted footprint drag friction load-intensity distribution for Space Shuttle 
nose-gear tire subjected to an inflation pressure of 300 psi. For model 1 : p stat j c = 0.6; Adynamic = 0.51 ; e n = 1 .0E+12; 
= 2.0E+03; = 0.5. 
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(b) Lateral friction load intensity, normal load = 2000 lbf. 



Footprint length, in. 



-4 


i 1 r 

-2 0 


2 


Footprint length, in. 


1 

4 


Load 

intensity, psi 


35 

25 

15 

5 

-5 


-15 


-25 


-35 



H 20 


10 


0 


-10 

-20 


-30 


-40 


(c) Drag friction load intensity, normal load = 2000 lbf. 


Figure 29. Effect of contact-area grid refinement on predicted contact load-intensity distributions. Space Shuttle 
nose-gear tire subjected to an inflation pressure of 300 psi and unsymmetric static loading conditions. For model 1: 
Fstatic — 0.6, Fdynamic ~ 0-51; £„ = 1.0E+12; E ; = 2.0E+03; E rc j ax = 0.5. 
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(c) Drag friction load intensity, normal load = 2000 lbf. 

Figure 30. Comparison of circumferentially-grooved tread tire and smooth-tread tire contact load intensities. Space 
Shuttle nose-gear tire subjected to an inflation pressure of 300 psi and symmetric static loading conditions. For 
model 1 : ^ static = 0.6; Hd ynam i c = 0.51 ; e„ = 1 .0E+12; e, = 2.0E+03; e^ax = 0.5. 
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Figure 31. Comparison of normal load intensities from frictional and frictionless contact. Space Shuttle nose-gear 
tire subjected to an inflation pressure of 300 psi and symmetric static loading conditions. For model 1 : p static = 0.6; 
^dynamic — 0-51, £ n = 1.0E+12; £, = 2.0E+03; £rg] ax = 0.5. 
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(a) Total strain energy density. 
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(b) Transverse-shear strain energy density. 


Figure 32. Variation of calculated strain energy density distribution. Space Shuttle nose-gear tire subjected to an 
inflation pressure of 300 psi and symmetric static loading conditions. For model 1: p static = 0.6; Adynamic = 0.51; 
e n = 1.0E+12; e, = 2.0E+03; = 0.5. 
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